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Abstract 

We discuss the possibiUty of measuring leading neutron production at the LHC. 
These data could be used to extract from it n^p and tt+tt^ cross-sections. In this 
note we give some estimates for the case of elastic cross-sections and discuss related 
problems and prospects. 
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1 Introduction 

In a recent paper [1] we pushed forward the idea of using the Zero Degree Calorimeters, 
ZDCs, designed for different uses at several of the LHC collaborations, to extract the total 
cross-section of the processes Trp — )• X and tttt — ;■ X at energies about 1-5 TeV. This would 
allow us to use the LHC as a vrp and tttt collider. In Ref. [1] there was also mentioned 
possible measurements of the elastic irp and tttt scattering. The physics motivation for 
extending this program to elastic scattering is very clear since the total and elastic cross- 
sections are so tightly interrelated (e.g., via unitarity) that any testing of various models 
of high-energy hadron interactions is little informative use without both of them (Fig. [1]). 
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Figure 1: a) Total pp cross-sections in the energy range 0.5 TeV < y/s < 15 TeV for 
parametrizations f^ (solid),l3^ (dashed), J^ (dotted) and IEI,IEI (dash- dotted). Data points are 
taken from Q/; b) evolution of the diffractive pattern in the pp scattering for ^/s = 1 — )• 20 TeV 
(from right to left) for parametrizations from J^ (solid) and ^8] (dashed). 

We would like also to remind to the reader an idea of universality of strong interactions 
at superhigh energies, in the sense that any ratio of two total or elastic cross-sections will 
asymptotically approach 1 independently of initial states (Fig. [2^). How close we are 
to reaching the asymptotic regime, asymptopia, can be tested when looking at the ratio 
of elastic to total cross-sections. Most of theories predict it to be 1/2 [2]- [6] (Fig. [Jt) 
though there are some theories which predict it to be equal to 1 [9], [10] (Fig. [St). Sure, 
this can and will be done. The TOTEM experiment has been designed to measure this 
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Figure 2: a) Ratio ofir^p to pp total cross-sections for COMPETE IS] parametrization; b) ratio 
of elastic to total cross-sections for pp, ir^p and vr^vr^ processes for parametrizations ^ (solid) 
and IEI,IEI (dashed); c) ratio of proton-proton elastic to total cross-sections in f^,|-?0|/. 
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Figure 3: Evolution of dijfractive pattern for pp (solid), Tr~^p (dashed) andTT~^TT~^ (dotted) elastic 
processes at a) \/s = 5 TeV and h) ^/s = 14 TeV for the parametrization from the Appendix B. 

ratio but any information about such a feature for the interaction of the lightest hadrons 
is impossible to overestimate. 

There are many other questions to be asked. For instance, how different are the interaction 
radii in pp,7Tp and vrvr high-energy colhsions? The properties of the interaction region could 
be obtained from diffractive patterns, which are different for these processes at the same 
energy (Fig. [3]). It would be interesting to know the dependence of these interactions on 
the very different quark-gluon content of colliding particles etc. If it will be possible to 
provide such a marvellous opportunity as access to the data on Tcp and vrvr TeV-energy 
interactions all the landscape of the "soft physics" will be transformed to the better. In 
what follows we present our estimates on the possibility to extract vrp and vrvr elastic 
cross-sections using CMS as an example. Just before sending this article to arxiv a paper 
on the four-body reaction pp — )■ nn7r~^n~^ at the LHC was placed there. However, the 
authors of [H] do not estimate possibilities for extraction of the tttt elastic cross-section. 

2 Exclusive single pion exchange 

The diagram of the exclusive single-pion exchange (SttE) process p + p^n + n^+p is 
presented in Fig.Hli. The momenta are pi, p2, Pn, v'-k-, v'2 respectively. In the center-of-mass 
frame these can be represented as follows (arrows denote transverse momenta): 
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Figure 4: Diagrams for the exclusive signal and background processes in leading neutrons pro- 
duction (initial rescattering corrections are not shown), a) Signal for the elastic ir'^p scattering: 
process with an exclusive single pion exchange (SirE) p + p — )• n + vr"*" + p, M is the mass of 
the vr+p system; h) background for the elastic n^p scattering: low mass single dissociation with 
Pomeron and reggeon exchanges; c) signal for the elastic 7r"'"7r"'' scattering: process with an ex- 
clusive double pion exchange (DttE) p + p — )• n + vr"*" + tt"^ -\- n, M is the mass of the tt+tt"'" 
system; d) background for the elastic tt'^tt^ scattering: SttE with single low mass dissociation in 
the TT^p channel; e) background for the elastic vr+vr^ scattering: low mass double dissociation 
with Pomeron and reggeon exchanges. 
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Figure 5: Amplitudes squared and cross-sections of the processes: a)p + p^n + X (SttE), b) 
p + p^n + X + n (DttE). S represents soft rescattering corrections. In this note X = tt^ + p 
(elastic SttE) and X = 7r+ + 7r+ (elastic DttE). 



As an approximation for tt exchange we use the formula shown graphically in Fig. [5^) . If 
we take into account absorptive corrections this formula can be rewritten as 
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where $x is the phase space for the system X produced in the 7i~^p scattering, the pion 



trajectory is aT^it) = a'^{t 



m„ 



The slope a' ~ 0.9 GeV ^, ,^ = 1 — xl, were xl is 



the fraction of the initial proton's longitudinal momentum carried by the neutron, and 
G^0p„/(47r) = C^+pn/lSvr) = 13.75 [12]. The form factor F(t) is usually expressed as an 
exponential 

F{t)=exp{bt), (8) 

where, from recent data [13], [11], we expect b ~ 0.3 GeV^^. We are interested in the 
kinematical range 0.01 GeV^ < |t| < 0.5 GeV^, ^ < 0.4, where formula ([7]) dominates 
according to [15j and [16]. At high energies we can use any adequate parametrizations 
of different vr+p cross-sections. Here we replace dax,-K+p/d^x in © by dad^Tr+p/dtp or 
integrated crei,n+p instead of crtot,7r+p in [1] • 

The suppression factor S arises from absorptive corrections [17]. We estimate ab- 
sorption in the initial state for inclusive reactions and for both initial and final states in 
exclusive exchanges. For this task we use our model with 3 Pomeron trajectories [8]: 

ajp^{t)-l = (0.0578 ±0.002) + (0.5596 ±0.0078)t, 

aip^ (t) - 1 = (0.1669 ± 0.0012) + (0.2733 ± 0.0056)t , 

ajp,{t) - 1 = (0.2032 ± 0.0041) + (0.0937 ± 0.0029)t , (9) 

These trajectories are the result of a 20 parameter fit of the total and differential cross- 
sections in the region 

0.01 GeV^ < \t\ < 14 GeV^ 8 GeV < y/s < 1800 GeV. 

Although the x^/d-o.f. = 2.74 is rather large, the model gives good predictions for the 
elastic scattering (especially in the low-t region with y^ /d.o.f. ~ 1). 

We use the procedure described in [18],[19] to estimate the absorptive corrections. 
With an effective factorized form of the, here-under, expression ( ITOl) used for convenience, 
we obtain: 

— ^^^^^^^ — = Ke +g )iMe,g )i (^3^ — ^^ — 5(./.o,e,g ), (10) 

(m2e + g2)i$5(e,g2)|2 • ^ ^ 

The functions $o and $s arise from different spin contributions to the amplitude 

Ap^n = -J===^n (mp^ ^3 ■ $0 + g <? ■ $s) ^p (12) 

and both are equal to $b in the Born approximation. Here (Tj are Pauli matrices and 
^n, ^p are neutron and proton spinors. Functions $o,s,b are given in the Appendix A. 
For the 7r~^p elastic cross-sections we use parametrizations from [4] and |5j,[6j which are 
described in Appendices B and C correspondingly. 

The differential cross-sections for the process p + p^n + n'^+p a.t ^/s = 10 TeV are 
depicted in Fig. [61 The total cross sections are listed in Table [H They are in the range 
8-270 fih for all values of ^max < 0.4 implying that we will have plenty of rate for the 
measurements. 



Table 1: Total p + p ^ n + tt'^ + p cross-sections in the kinematical region < |^ < 0.5 GeV, 
^min = 10~^ < ^ < (,max for two parametrizations given in Appendices B (C). 



<,max 


0.05 


0.1 


0.2 


0.3 


0.4 


C'p+p— i>n+7r++p ; A'-t) 


8.5 (8.8) 


37.5 (39) 


128 (132) 


208 (214) 


259 (266) 




b) 



dff 


mb 


dt' GeV^ 




5- 


1.5 


/>*. 


1.25 


/ \ 


1 


/ \ 


0.75 


/ \ 


0.5 


/ ^V 


0.25 


' ^--^ 



0.05 0.1 0.15 0.2 0.25 



Figure 6: Integrated cross-sections of the SnE process p + p — ;• n + vr"*" -\-p for parametrizations 
from Appendices B (solid) and C (dashed), a) ^ (13.01 GeV^ < |t| < 0.5 GeV^;,- b) ^, lO^^ < 
e < 0.3. 

At low energies {^/s < 70 GeV) the region of applicability of our model is given by the 
inequalities 

0.01 GeV^ < |t| < 0.5 GeV^ 10"^ < ^ < 0.4. (13) 

At higher energies this region may be smaller (say ^ < 0.2), since this corresponds to 
masses M = 3 TeV at ^/s = 10 TeV, and for larger masses the approach may break down. 



3 Exclusive double pion exchange 

As noted above, the Double pion Exchange (DttE) process can give information on both 
total and elastic vrvr sross-sections. vrvr cross-sections have been extracted in the past using 
the exclusive cross-section [2U] . The results are shown in Fig. [71 There is some tendency 
for an early flattening of the tttt cross-sections. In irp and pp cross-sections this flattening 
begins at higher energies and precedes further growth. 

We can extend the analysis for one pion exchange described above to double pion 
exchange (Fig. [5b, DvrE). The kinematics of the exclusive DvrE {p + p -^ n + n^ + ii^ + n, 
the momenta are pi, p2, Pm, p'^y P^2' P"-^ respectively) is similar to the exclusive double 
pomeron exchange process: 



Ptt. 



^.{-ir^^A 



Pn, =Pi- Pn, 



(14) 



M'^ = (p^i + P^a) - ^i6s - (gi + q2) ^ 66s, 



-t 



2 _ /„ , ^ \2 
7r2y 



'-TTTT KPtTI Ptti 



I \2 



(15) 
(16) 



n n cross-sectjotts in mb 



31- p — * 71- n+ n 



wo 

HI 

9 



■ n tutiil 
x G elastic 



( h 



n 



*'tttHllt 



ft 



SO: 
9 



luiuii L 



l.a 2.D S.O 4.0 

MCn-n+), GeV 



n- p — > 31- ;i- A-H- 



Mm 



M 






■Mlii 



lUU 



1,0 2.0 3.0 4,0 
M(]i-Jt-), GeV 



Figure 7: Elastic and total cross-sections for vr tt^ anc? vr vr scattering from the data on 
exclusive reactions as a function of the dipion invariant mass (Fig. 5 from Ref. !i20f). 



The cross-section can be evaluated as follows: 
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For the ti^tt^ elastic scattering we get 
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where functions p^j, ^jj and $b are given in the Appendix A. 

We are now ready to make predictions for high energies. Numerically calculated cross- 
sections for the exclusive DttE are shown in Figs. [H] and listed in Table |2] for parametriza- 
tions from Appendices B and C 



Table 2: Total p+p — )■ n+7r"'"+7r++n cross-sections in the kinematical region < |g| < 0.5 GeV, 
Cmin = 10^"^ < ^ < &nax for two paramctrizations from Appendices B (C). 
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Figure 8: Partially integrated cross-sections of the DttE process p + p — )• n + vr"'" + vr"*" + n at 
y/s = 10 TeV for parametrizations from Appendices B (solid) and C (dashed): a) dcr/d^id^2 for 



^1 = .^2 = ■? and < 1^1,2! < 0.5 GeV; b) da/dqidq-/ for \qi 



ml 



\^ and W-^ < 6,2 < 0.3. 



4 Extraction oi n^p and n^n^ cross-sections. The role 
of absorption. Backgrounds. 

To extract n'^p and 7r+7r+ cross-sections from the SttE and DvrE processes we can use 
equations ([7]) and ( 1T7|) . Let us rewrite these two equations in the following form 
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An exact extraction procedure is quite delicate. If we want to extract vr+p and tt+tt^ 
cross-sections in a model independent way, we have to take equations f l20|) . fl2Tl) in the 
limit ti — )■ 171^. For this limit we should extrapolate the parametrizations of the data 
on SttE and DttE differential cross-sections to the positive ti = rri^, i.e. beyond the 
physical region. Functions S and 5*2 are equal to unity for this value of ti, that is why the 
phenomeno logical model for these functions is not important. This procedure is actually 
the Chew-Low extrapolation method j21|.|22j. 

Experimentally extrapolation to m^ is rather difficult (see section[5]), since the errors in 
t are larger than m^. To get around this problem we extract cross-sections for pions with 
low virtualities and assume that the values (12(I]) .( pT]) are close to reality. It is clear from the 
fact that the main contribution to the cross-section comes from the region |tj| < 0.25 GeV^ 
(see Fig. |6]d). In this region the dependence of a^^+p {o'tt+tt+) on t is assumed to be weak 
enough. 
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Figure 9: Function 5(s/so, ^, qt) at y/s = 10 TeV in the physical region of negative t values for 
three different fixed (^ values: ^ = 0.3 (dotted), ^ = 0.1 (dashed) and ^ = 10~^ (solid). For low 
^ and \q\ function S is close to unity. 
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Figure 10: Function 52(s/so,^i,2, |9i,2|) in the physical region of negative t values at ^/s = 
10 TeV for: a) for fixed 6,2 = 0.01 h) for fixed |gi,2| ~ 0. c) Function F(6,6) aty/s = lQ TeV. 



Functions 5* and 5*2 are close to unity in the physical region of negative t values (see 
Figs. |9|T0|) . and we can estimate errors of the model due to absorptive corrections. It was 
shown in pLj that such a model dependent extraction works satisfactory for i/i < 70 GeV. 
This fact is illustrated in Fig. [TTb . All the parametrizations are close to the extracted 
values and the real data points, but for higher energies (Fig. [TTb) the difference between 
models becomes larger. 

To avoid singularities in the extrapolation procedure at t = and model dependence in 
5* and 5*2 it is convenient to extrapolate quantities in the r.h.s of fl20|) and fl2T|) multiplied 
by S t/m^ and 5*2 tit2/m^ correspondingly. The behaviour of S t/ml. is shown in Fig. [T2l 
It is a smooth function of t in the whole region of the extrapolation. Practically we will 



mz 



have c^TT+pS t/ml as a result of the extrapolation, which is equal to cr-^+p at t 

The role of absorptive effects (i.e. model dependence of the final result) is significant if we 
want to extract np and vrvr cross-sections from the SvrE and DvrE differential cross-sections 
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Figure 11: Total ■n'^p cross-sections versus different parametrizations: f^ (solid), fSi^ (dashed),^ 
(dotted) and f5^,f^ (dash- dotted), a) real data from PDG (triangles) up to ^/s = 25 GeV and 
extracted values (boxes) up to y/s = 70 GeV (see fB/); b) total -ir~^p cross-sections in the energy 
range 0.5 TeV < a/s < 7 TeV. 
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Figure 12: Function S{£,,t) t/rn^ versus t/m'^ at fixed ^ = 0.05. The boundary of the physical 
region to = — fn^^^/(l — ^) is represented by vertical dashed line in b). 

integrated in the wide region of t values, where absorption is strong. That is why we need 
an experimental instrument to measure differential cross-sections for low t values with 
good resolution. The present design of detectors does not allow t measuremets, it gives 
only restriction |t| <~ 1.2 GeV^ \i\. If to assume a weak enough t-dependence of Tip and 
TTTT cross-sections, then we could hope to extract these cross-sections (though, with big 
errors) by the following procedure: 
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Figure 13: Values of absorptive corrections integrated with form-factors in the region 
0.01 GeV^ < |t,| < 1.2 GeV^. a) S{^); b) 52(6,6)- 6 = 6 (solid), 6 = 0.1 (dashed), 
6 = 0.2 (dotted) and 6 = 0.3 (dash- dotted). 

Table 3: Total absorptive corrections for exclusive SnE and DttE in the kinematical region 



0.01 GeV^ < \ti\ < 0.5 GeV% 6 
B,C. 



10 < 6 < Cmax for parametrizations from Appendices 



<,max 


0.05 


0.1 


0.2 


0.3 


0.4 


•Jtot, SttE 


0.39 


0.34 


0.27 


0.22 


0.18 


•Jtot, DttE 


0.47 


0.41 


0.33 


0.26 


0.25 



Functions 5*2(656) and 5(6 are depicted in Fig. [I3l To suppress theoretical errors of S 
and ^2 we have to measure total and elastic pp rates at 10 TeV, since all the models for 
absorptive corrections are normalized to pp cross-sections. At present we can estimate 
the theoretical error to be less than 10% for this method from predicted values of total 
pp cross-sections in the most popular models (see Fig. [T^). 

The case of the DttE is more complicated since the function ^2 is not factorisable. For 
low ti it is approximately equal to 

F(6,6) = ^2(5/^0,6,6,0,0)^ 

V5(s/so, 6,0) + V5(s/so, 6, 0) - V5(s/so, 6, 0)5(s/so, 6, 0)) ' , (25) 



which is clear from Figs. ITOb.c. 

To estimate total absorptive effect we have to take the ratio 




n': e 



lO"'' < 6 < 6iax, 0.01 GeV^ < |i| < 0.5 GeV" 



(26) 
(27) 



$ is the phase space for the SttE (DttE), and dao/d^ is the cross-section without absorptive 
corrections (i.e. for S = S2 = I). The results are listed in Table [3l 

If we could measure momenta of all the final particles, we would have only exclusive 
backgrounds for our SvrE and DttE signal processes with elastic n^p (Fig. H^) and 7r+7r+ 
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(Fig. Ht) scattering. The single low mass dissociative background for the exclusive SvrE 
is depicted in the Fig. |3)d. For the exclusive DvrE we have two exclusive backgrounds 
of Figs. |Hi,e. In a real experiment we can detect only one or two particles in the final 
state, and we have to take into account all the inclusive backgrounds: single and double 
dissociation, central diffraction, minimum bias with neutrons production, inclusive single 
and double charge exchanges with 7r+ (see Fig. [T^ and also p, 02 exchanges (for more 
exact estimations of inclusive backgrounds see Section E]). 

5 Experimental possibilities 



In this chapter we analyse CMS [23] capabilities to measure elastic vr+p and 7r+7r+ scat- 
tering at 10 TeV, cm. energy of LHC protons in the first runs. The CMS Zero De- 
gree Calorimeters ,ZDCs, [2U [25] can measure leading neutrons in the exclusive SttE, 
pp — )■ npi:^ (Fig- SK), and DvrE, pp — )• nvr+vr'^n (Fig. HJ:), processeqj. The ZDCs are 
located between the two beam pipes at 140 m on each side of the interaction point. 
They are able to measure the energy of forward neutral particles in the pseudorapidity 
region \ri\ > 8.5. SvrE and DvrE events have been generated in the framework of the 
simulation package EDDE [26] . The kinematics of the SvrE and DvrE proceses are defined 
by ^„ and t„ of the leading neutron. The Ptt^^j^ vertex is generated on the basis of the 
model described in Ref. [^. For the simulation of tt^^P ^^^ '^virt'^virt elastic scattering 
PYTHIA 6.420 [27] has been used. Obviously, vrj^^^ and p (vrj^^ ^^^ vrj^^i) ^^^ interact 
inelasticaly and diffract ively. Then, in the diffractive interaction of ttJ^^^ and p, either the 
''^virt or the p, or both of them, can dissociate. All of these processes have been studied 
as backgrounds, as well as minimum bias and diffractive pp events. Diagrams for some 
of the background processes are shown in the Fig. [TH Signal and background have been 



^ Further, exclusive SttE (DttE) elastic events, i.e. SttE (DttE) with TTyi^tP i'^tirt'^virt) scattering 
elasticaly, will be designated as Sir'Eeiastic (D'l'Eeiostic) for brevity. 



y^ X X ^ 

b) 





Figure 14: Inclusive SttE (a,h,c) and DttE (d,e,f) backgrounds, a) single dissociation, h) dou- 
ble dissociation, c) minimum, bias in the 7r~^p channel, and d) single dissociation, e) double 
dissociation, f) minimum bias in the 7r"'"7r''' channel. 
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generated by EDDE v. 3. 0.0 and PYTHIA 6.420. The SvrE cross section, including all 
types of Ti^irtP interactions, is estimated to be about 2.6 mb at 10 TeV and ^„ < 0.4 [T]. 
Corresponding cross sections for the signal and backgrounds are listed below. 



Signal: 



SvrEeiastic '■ <7j; 



0.33 mb. 



Backgrounds from pp and inelastic SvrE events: 

- minimum bias events: <7pp^x = 50 mb; 

- single diffractive dissociation: app^px = 14 mb ; 

- double diffractive dissociation: app^xY = 9.7 mb; 

- SvrE, minimum bias in the vrJ^^iP channel. Fig. [Mb : app-^nX = 1-54 mb; 

- SvrE, single diffraction in the Ti^i^tP channel with proton dissociation. Fig. [T^ : 
o'pp^mr+x = 0.23 mb; 

- SvrE, single diffraction in the Trj^^tP channel with 7r+ dissociation. Fig. [T^ : 
CTpp-^ripX = 0.20 mb; 

- SvrE, double diffraction in the vrJ^^tP channel. Fig. [T4b : CTpp^nXY = 0.27 mb. 

Fig. [m shows the ratio of events for SirEeiasUc (shadowed) and background processes. 
On the picture minimum bias processes have numbers less than 90 and diffractive processes 
have numbers 92, 93 and 94, according to PYTHIA's definition. All SvrE processes are 
placed between numbers 400 and 500 and DttE processes are in the region from 500 to 
600 in our generation. Signal processes, SnEeiasUc and DnEeiasUc, have numbers 491 and 
591 respectively. SnEeiasUc events contribute ~0.4% to the total cross sections, DvrEeZastic 
events is around 0.025% only. 

DttE events have been simulated by the same method, as SvrE, using EDDE v.3.0.0 
and PYTHIA 6.420. DttE cross section, including all types of vrj^^^irt interactions, is 
estimated to be about 200 /ib at 10 TeV and ^ni,n2 < 0.4 [1]. Corresponding cross sections 
for the signal and inelastic DttE backgrounds are listed below. 




30 \ 200 

DilTruclian 



300 
N of process 



Figure 15: Events ratio for S7rE(,iastic (shadowed) and background. 
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Signal: 

Backgrounds from inelastic DttE events: 

- DvrE, minimum bias in the T^tirt'^tirt channel, Fig. IT^f): (Jpp^nXn = 124 /ib; 

- DvrE, single diffraction in the vr^^^7rj".j,j channel, Fig. [HJi: app^n^+xn = 30 /ib; 

- DttE, double diffraction in the vrj^rt^iliri channel. Fig. [T4re): (Jpp^nXYn = 22 /ib. 

Inelastic and diffractive pp interactions produce background events for DirEeiasUc, as well 
as for SnEeiastic- Moreover, SvrE elastic and inelastic processes produce strong back- 
grounds for DnEeiastic, in addition to the inelastic pp and DvrE. And, on the contrary, 
DttE elastic and inelastic events can imitate SnEeiasUc- In our simulation, cross sections 
for the SirEeiastic and DnEeiasUc signals, depend on vrj^^jP and vr^^'^'^lirt elastic scattering 
models integrated to the PYTHIA 6.420. It is interesting to note that values obtained for 
these cross sections are very close to those ones which can be calculated in the BSW [4J 
and GP [HI E] parametrizations, see Table [1] and [2l The ratio between the SttE elastic 
(signal) and inelastic (part of background) events are presented in Fig. [T6k . Fig. [T6b 
shows the same for DttE. 

As in paper [1], for SvrE selections we choose events with neutrons in the forward or 
backward ZDC and with the absence of neutrons in the opposite one: 



N^ > & N^ = 
N^ > & Nl = 0. 



(28) 



For the DvrE, we selected events with neutrons in both the forward and backward ZDCs: 



K>o 



& 



N^>0. 



(29) 



Here, A^^ (A^n) is the number of neutrons hitting the forward (backward) ZDC. Such 
selections suppress the background for SnEeiasUc (^''^^eiasUc) events by a factor 14 (160). 



o) 



b) 



> 

a; 
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DttE ioi«i 



. D^rKelaslic 




WTeV 



Figure 16: (a) SirE elastic (solid) and total (dotted) events distribution versus the (tt^^j-^p) in- 
variant mass; (b) DttE elastic (solid) and total (dotted) events distribution versus the ('J^^irt'^tirt) 
invariant mass. 
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The signal to background ratio becomes equal to 0.05 for the SnEeiasUc and 0.04 for the 
DTiEeiastic- Fig. [T7] shows the ratio of events for SnEeiasUc (shadowed) and background 
processes after the selection (1251) . The same picture is plotted for the DnEeiasUc signal 
and background after the selection f l29|) . Fig. [181 

The signal SnEeiasUc {pp —> mc'^p) event has neutron, proton and vr"*" in the final state. 
Apart from neutrons, which can be detected by ZDC, two other particles move out of the 
CMS acceptance. The proton, scattered elasticaly, should move inside the beam pipe. 
The 7r+ meson should fly in the same direction as the neutron and it is deflected at a 
small angle too (see Fig. I20b). Thus, though it looks as a paradox, we should demand an 
absence of a signal in the CMS detectors, except of the one of ZDCs, for the STrEeiasUc 
trigger. For example, we could select events with zero signal in the CMS calorimeters: 



Nbarrel 
Nendcap 
Nhf = 

NcASTOR 






0. 



(30) 



Fig. Hn] shows the efficiency of such selection for the SttE elastic- The signal to background 
ratio becomes equal to 1, i.e. we have achieved improvement 18 times better in comparison 
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Figure 17: The ratio of events for the signal SnEeiasUc (shadowed) and background processes 
after the selection 
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Figure 18: The ratio of events for the signal DirEeiasUc (shadowed) and background processes 
after the selection 
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Figure 19: The ratio of events for the signal SirEeiasUc (shadowed) and background processes 
after selections ^ & (3^. 



O) 



1QQ0 r L| 
500 



KOJOeV 



I t I f-^ --' 



0,9 1 

U Gev' 



400 

200 





bl 



_ 




- f. 

■* y I. 


n>9 


I-I 

■P 

ill' 





10 



^QQ 



>00 



7.5 



Figure 20: a) t of the leading neutron for the signal STTE(,iasUc (solid) and background (dotted) 
after selections i28\) & i30\} : b) r/ of the 7r+ for the signal SirEeiastic (solid) and background 
(dotted) after selections [28\) &: (3^ k. tn < 0.2 GeV^; c) the (ir^p) mass distribution for the 
signal S-KEeiasUc (solid) and background (dotted) after selections ^28\) k, (W) k. tn < 0.2 GeV 
& ri^+ > 8.5. 



with the previous selection f l28|) . Detailed study of the rest of backgrounds (left bin on 
the Fig. [T9|) has shown that it contains processes of the single diffractive dissociation, 
pp — )■ pN* — > priTc^, where the Pomeron exchange leads to the proton excitation A^* and 
its subsequent decay to the 7r+ meson and neutron (see the diagram on the Fig. [T^ . This 
reaction can imitate the SiiEeiasUc process as well. However, the further careful study of 
both reactions has shown some difference in their kinematics. Thus, t distributions of 
neutrons have different slope parameters. It could improve the signal/background ratio 
up to the value 1.7 by the selection |t„| < 0.2 GeV^, Fig. [20k . 

Each of the CMS ZDCs consist of two sections, an electromagnetic, EM, part for mea- 
suring photons, ttq, rjs etc and a hadronic part designed to measure neutral hadrons such as 
neutrons and As [25]. The energy resolution of the detector for hadrons is 138%/\/E+13%. 
The electromagnetic part is divided into strips that run in the vertical direction. These 
strips can be used to measure the horizontal position of the particle's impact point with 
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Figure 21: dN/d^dr for a) ^s = 900 GeV and b) 7 TeV 
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Figure 22: The t value for neutrons as a function of distance from the coUision axis and 
i for a) yi = 900 GeV and b) 7 TeV 



a resolution of about 0.4cni. The hadronic part that has is divided into 4 depth seg- 
ments but has no transverse segmentation. About 1/3^'°' of the time neutrons will start to 
shower in the electromagnetic part and for these neutrons we can extract some position 
information. The ZDCs can also be used to select events in the CMS level one trigger. 

The geometrical acceptance of the calorimeter is ±4.4 cm horizontally and ±5.0 verti- 
cally. For CMS the LHC beams cross in the horizontal plane and so the nominal position 
of the zero degree point will vary depending on the crossing angle. For example if the 
crossing angle of the beam is 140/iradians the zero degree point will be at x=±2cm. Given 
the energy and position resolution of the detector it may be possible to make a rough 
measurement of the angular distribution of the neutrons. 

The only independent measurements that the ZDCs can make of the neutrons are the 
energy loss ^ and the distance from the collision axis r = a/x^ ± y^. Figures |2T1 122] and 
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Vs'=0.9 TeV 



VF^T TeV 




a) 



b) 




Figure 23: The relative resolution for t, i.e. 6t/t value for neutrons as a function of 
distance from the collision axis and ^ for a) ^/s = 900 GeV and b) 7 TeV. We have 
assumed that 6r = 0.5cm, distance from the detector to the interaction point is 140 m 
and 6^/^ ~ 0.14 



show the distribution of neutrons, their t value and the resolution of t versus ^ and 
r for a) ^/s = 900 GeV and b) 7 TeV. As the energy increases the radial distribution 
tends to shrink towards r = 0. Given the current limited position resolution of the ZDC 
it may be possible to gain some information about the t distribution at ^/s = 900 GeV. 
However for multi TeV energies it will probably be necessary to upgrade the detector. 
Figure EOj^b) shows that the vr^ mesons are deflected a bit stronger for the signal than for 
a background. If we could have counters for charge particles in the pseudorapidity region 
T] > 9, it would allow us to improve the signal/background ratio up to 5 and higher. 
Fig. I20b . The possibility of such counters installation along the LHC beam on both sides 
of the CMS was studied in the Ref. [28]. Set of FSCs, placed at distances from 60 to 
140 m from the interaction point, could cover the pseudorapidity region from 8 to 11. 
They could register particle showers induced by the primary tt"*" with high efiicency, up 
to 70%. Unfortunately, in the present setup of the forward CMS detectors there are no 
forward counters and present design of the ZDC does not allow to measure t of the leading 
neutron. So, this is a task for the future. 

Selections fl30|) . applied to the DnEeiasUc, improve the signal/background ratio up to 
1.1 (see Fig. 1241) . The rest of the background comes from the the double diffractive 
dissociation pp — )■ N*N* — > mc^Tc^n (the left diagram on the Fig. ^M and from the SvrE 
events produced by single diffraction in the ir'^p channel with the subsequent decay of the 
excited protons to vr"*" and neutron (the right diagram on the Fig. [2^ . t distributions for 
signal and background are different, as it is shown on the Fig. I25k . As for the SvrE, we 
could improve the signal/background ratio up to ~ 2 by selections: 



\4\ < 0.25 GeV^ 
\t\ <0.25GeV2 



(31) 
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Figure 24: The ratio of events for the signal DTiEf.iastic (shadowed) and background processes 
after selections (Eg) & IM) . 
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Figure 25: a j t of the leading neutrons for the signal DnEf^iasUc (solid) and background (dotted) 
after selections ( fIP|) & f3U\): b) rj of the vr^ mesons for the signal DirEeiasUc (solid) and back- 
ground (dotted) after selections i flPj) &: [3U\) &: |g7]) ; c) the (V+vr+j mass distribution for the 
signal DnEeiasUc (solid) and background (dotted) after selections [2^) &: (W) & [3l\) & (W) . 



using f l3T|) in a combination with fl29|) and fl30|) . For further improvements of data we 
could use different deviation of tt meson from the beam axis for signal and background, 
Fig. [25b . Selections 

give a data sample with a signal/background ratio ~ 7, see Fig. 



(32) 



6 Discussions and conclusions 

In conclusion, our study of '^n'Ei elastic and DvrEe/astJc processes shows that with present 
setup of the forward detectors we could expect observation of STrEg^astic and 'Dti'Ej elastic 
events mixing with background in the proportion ~1:1. Rough estimations on the gen- 
erator level shows that we could observe ~ 10^ STrEe^astic events distributed in the mass 
region from 1 to 6 TeV and ~ 10^ 'DTi'EeiasUc events distributed in the mass region from 
0.5 to 4 TeV at the integrated luminosity 1 pb~^. As it was said they will be mixed with 
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approximately the same amount of background events. Improvement of the data purity 
demands a considerable modernisation of the forward detectors. Some modification of 
the ZDC is required to measure the t of the leading neutrons. It would be very useful 
to install forward shower counters FSCs along the beam at distancies from 60 to 140 m 
for the detection of elastic scattered vr-mesons in the region ?] > 8 would improve the 
measurements significantly. Realisation of such modifications is beyond this article. 

Theoretically, it is very interesting to have both elastic and total cross-sections of Tip 
and TTTT scattering. At present we could only use the extraction procedure for t-integrated 
SttE and DvrE cross-sections ( 123|) . (!24|) which is far from the ideal one ( 120|) . (pTj) . Rough 
estimations give the model error about 10%. The main part of this error comes from the 
uncertainties in the absorptive corrections which are normalized to pp total and elastic 
cross sections. Measurements of the pp total cross section, which would be done by the 
TOTEM experiment at LHC, can improve the precision of our model-dependent extraction 
procedure significantly. We would like to stress again, that model-independent extraction 
procedure for np and tttt total and elastic cross sections makes precision measurements of 
t of the leading neutron at small angles mandatory. 

In spite of all the difficulties, proposed tasks are of exceptional importance, and we 
hope that they will push modernisation of the forward detectors for future precise mea- 
surements. 

Appendix A 

Here we define functions for the calculation of absorptive corrections for the SvrE and the 
DttE processes. For the SvrE we have: 

- ^^-^T-k^' ^"^ 0' (33) 



27r g2 + e2l + /32g; 
iV(0 = (l-0%^e^exp' 



1-e 



(34) 



6 + a' In i 
/3^ = ^^' '' = <^" + ^' (1 - 0, (35) 

oo 

% = ^JdbQo{b,^Amib), (36) 





oo 



ms = ^JdbQs{b,^,mv{b), (37) 





bJo{b\q\)(Ko{eb)-Ko(^, 
0o(6,e, m = T^^^ ^^^' (38) 



6Ji(6|g1)(eiri(e6)-ii^i(|)) 

1 - /32e2 



0s(&, e, m = \ ,, , ^^^, (39) 
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Table 4: Parameters of the model I8j. 



i 


1 


2 


3 


Ci 


53.0 ±0.8 


9.68 ±0.16 


1.67 ±0.07 


rf (GeV-2) 


6.3096 ±0.2522 


3.1097 ±0.1817 


2.4771 ± 0.0964 



V{b) = exp (-fiei(s/so, b)) , 

^ ^ 2c- / s 



IGnBi V So 



a/p,(0)-l 



exp 



4a 



•^/P, 



Ei = a'_ In -e-*2 + 4 



So 



rp^ 



the values of parameters can be found in ([9]) and in Table HI 
For the DttE process we can write the following expressions: 



(40) 

(41) 

(42) 



^^^ = ^^^^^^j^5'^ j dhdhQ,{h,i,, |gi|)e,(&2,6, \q2\)h{hM), 



TT 

hiPl. &2) = J'^y Ubl + bl-2b,b2COsA , 


Poo = mlC,iC,2, Pos = rripC,!, Pso = mp^2, Pss = 1- 



(43) 

(44) 
(45) 



Appendix B 

For the calculation of elastic cross-sections we use the Bourrely-Soffer-Wu (BSW) parametriza- 
tion |1]. Functions and values of parameters are given below. 



00 






(46) 
(47) 



fio(s,&) = njp + J2^i 



n 



IP 



In' s 



1 + 



(1 + nfi)' 



FBswbioT s >mp, |t|, 



(48) 
(49) 
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Table 5: Parameters of the model ^. 



c 


c' 


mi 


7712 


mzn 


u 


Ott 


/ 


a 


0.167 


0.748 


0.577225 


1.719896 


0.7665 


4.2414 


2.3272 


6.970913 


1.858442 



Table 6: Parameters of the model U^ for secondary reggeons in pp scattering. 



i 


bJ 


A2 


P 


hi, GeV"" 








8.54 


c. 


-167.3293 


-24.2686 


124.91969 


a^{t) 


0.3229 + 0.7954t 


0.3566 + t 


0.3202 + t 



For the vr+p elastic scattering we have i = p in fHHj) and 



al-q^ 1 



i^B/n^W = / g ^g Jo(g6)/.- ^ ^ + 4)(1 + 4 





X ap(0)-l -4^ 



where values of parameters are listed in Table [51 
For the pp elastic scattering i = u, A2, p and 



00 
FZwib) = J q dq Mqb) f 



a^ — q^ 



«^ + ^^ (1 + 5)^(1 + 5) 



fi.(6) = qdq Joiqb) Cie~''^'^\l ± e-'^'^d-i^)) ( A 




«,:(-g') 



(50) 



(51) 



Bp = bp + a'JO) In — , bp = 4.2704, ap(t) = 0.3202 + t, C^ = 4.1624, (52) 
So 



(53) 
(54) 



Values of parameters are listed in Table El 

In this paper we take the following parametrization of the ir'^n'^ elastic scattering, which 

is based on the BSW |1] one (approximate expressions for s ^ Sq, \t\ ^ 1 GeV^) 



TTTT ^ 



«L + g' (1 + ;V)2 



Fsswib)^ J qdqJoiqb)!: 


/^ 1 2 1 



f^o - fi/p, /. 



TTTT /* 5 



/ ' aL al a^ 



(55) 



(56) 
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Appendix C 

Another parametrization for pp, vr+p and vr+vr"'' cross-sections is taken from [5], [6]. The 
scattering amphtude is represented in the usual eikonal form 

p2iS{s,b) _ 1 

T{s,b) = (57) 

(here T(s, b) is the amphtude of the elastic scattering in the impact parameter b space, 
s is the invariant mass squared of colliding particles and S{s,b) is the eikonal function). 
Amplitudes in the impact parameter space and momentum one are related thorough the 
Fourier- Bessel transforms 

f{s, b) = -^ r d{~t)MbV^t)f{s, t) , (58) 

167rs Jq 

f{s,t)=A7rs db\Ub./^t)fis,b). (59) 

Jo 

Eikonal function in the momentum space is 

5{s,t) = 5p{s,t) + 5f{s,t) = 

The parametrization for the pomeron residue is 

/3p(t) = Bpe^^\l + dit + d2t^ + d-it^ + dit^) , (61) 

which is approximately (at low values of di, d2, d^ d^) an exponential at low t values. 
Residues of secondary reggeons we set as exponentials: 

Pf{t) = Bfe''fK (62) 

Phenomenological parametrization for the "soft" pomeron trajectory is set to 

ap{t) = 1 + pi l-p2t f arctg(p3 - P2t) - -j . (63) 

Trajectories of secondary reggeons /2 and u are parametrized by functions 

/ 3 \ 1/2 

«R(t)= f — «s(V-t + CR)j ,R=/,a;, (64) 



where 

" A2 -•- A2 
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Table 7: Values of parameters of the model 1^,IEI for the pp scattering. 





Pomeron 




/2-reggeon 




w-reggeon 


Pi 


0.123 


c/ 


0.1 GeV^ 


Cu, 


0.9 GeV2 


P2 


1.58 GeV-2 










Pz 


0.15 










Bp 


43.5 


5f 


153 


B^ 


46 


bp 


2.4 GeV-2 


6, 


4.7 GeV-2 


b^ 


5.6 GeV-2 


d. 


0.43 GeV-2 










d. 


0.39 GeV-4 










d. 


0.051 GeV~6 










d/j^ 


0.035 GeV-s 










ap(0) 


1.123 


«/(0) 


0.78 


a.(0) 


0.64 


«p(0) 


0.28 GeV-2 


a[,{0) 


0.63 GeV-2 


«L(0) 


0.07 GeV-2 



Table 8: Values of parameters of the model J^jJBi for the n^p scattering. 



n^p ^ TT^p 


Bp 


26.7 


bp 


2.36 GeV-2 


di 


0.38 GeV-2 


d2 


0.30 GeV-^ 


ds 


-0.078 GeV-6 


d. 


0.04 GeV-8 


Bf 


67 


bf 


1.88 GeV-2 



is the one-loop analytic QCD running coupling [29], Uf = 3 is the number of flavours, 
A = A*^^) = 0.346 GeV [30j. Parameters cj , c^ > are rather small to spoil the asymptotic 
behaviour of secondary trajectories in the perturbative domain. 
Residues for tttt, np and pp scattering are assumed to be 



(3Vit) 



PY{t) 






(66) 



(67) 



Parameters of the model are listed in Tables ITfHl 
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